NASA Technical Paper 1548 


NASA 

TP 

1548 

c .1 


LOAN COPY: i' : : j ’.jriT-I T'* v 
AvF'Vc L j t..c. • r I i\’ ! . L. ' l;‘ i' £/'• 

KIRI LAND AfB, N . M. 


The Wave Structures of the Eady 
Model of Baroclinic Instability 



Jae Min Hyun and William W. Fowlis 


OCTOBER 1979 





TECH LIBRARY KAFB, NM 


TECH LIBRARY KAFB, NM 



NASA Technical Paper 1548 


The Wave Structures of the Eady 
Model of Baroclinic Instability 


Jae Min Hyun and William W. Fowlis 
George C. Marshall Space Flight Center 
Marshall Space Flight Center, Alabama 


NASA 

National Aeronautics 
and Space Administration 

Scientific and Technical 
Information Branch 


1979 




TABLE OF CONTENTS 


Page 


INTRODUCTION 1 

FORMULATION AND EADY’S SOLUTION 2 

RESULTS AND DISCUSSION 9 

REFERENCES 21 



LI ST OF ILLUSTRATIONS 


Figure 

Title 

Page 

1 # 

Propagation speed C versus scaled wavenumber K- , * . . 

. 9 

2. 

K 

Growth rate KC^ versus scaled wavenumber K 

. 10 

3. 

Amplitude and phase variations of the most unstable 
mode ( K - 1. 6061, C R = 0. 5, Cj - 0. 1929) 

. 11 

4. 

The x-z plot showing the interrelationships among the 
dynamical variables 

. 12 

5. 

Amplitude and phase variations of an unstable wave 
(K- 2.2, C - 0.5, C T =0.0961) 

X\ i 

. 14 

6. 

Amplitude and phase variations of an unstable wave 
(K- 2.35, C = 0.5, C = 0.0476) 

. 15 

7. 

Amplitude and phase variations of a marginally stable 
most (K- 2.3996, C R = 0.5, Cj = 0.0) 

. 17 

8. 

Amplitude and phase variations of a fast-propagating 

stable mode (K= 3.5, C = 0.7131, C=0.0) ........ 

K I 

. 18 

9. 

Amplitude and phase variations of a slow-propagating 
stable mode (K= 3.5, C R = 0.2869, C = 0.0) 

. 19 


iv 



THE WAVE STRUCTURES OF THE EADY MODEL OF 
BAROCLINIC INSTABILITY 

INTRODUCTION 


A number of baroclinic fluid models have been proposed to investigate 
the large-scale (synoptic scale) processes in geophysical fluid dynamical sys- 
tems [1,2]. It is now generally understood that the initial development of the 
large-scale travelling waves in the westerlies in temperate latitudes may be 
modelled by use of a linear stability analysis of a structured zonal current of a 
baroclinic fluid. The fundamental source of energy for the growing disturbances 
is recognized to be the latitudinally differential solar heating, which creates a 
temperature gradient in the north-south direction. 

Under the assumptions of geostrophic balance and hydrostatic approxima- 
tion, which are appropriate for scales of long waves in the middle latitudes, 
systematic derivation of the scaled governing hydrodynamic equations of motion 
yields a set of approximate forms of equations known as the quasi-geostrophic 
equations [3], 

Among several models, the well-known model of Eady [2], using the 
Boussinesq fluid approximation and neglecting the /3-effect, allows ready 
analytical solutions; it captures the essential features of the baroclinic waves. 

In Eady’s model; a (vertically) stably stratified zonal current is contained 
between two rigid, parallel, rotating, horizontal boundaries. Two principal 
results of the Eady model are that there is a short wavelength cutoff for the 
unstable wave region and that the short wavelength cutoff and the most unstable 
wavelength do not depend on the vertical wind shear. 

In this report, the solutions of the Eady eigenvalue equation are recast 
into a somewhat more manageable form, and the propagation speed and the 
growth rate of baroclinic waves of varying sizes pertinent to the Eady model 
are computed. Another purpose of this report is to provide a comprehensive 
picture of the detailed wave structures, characterized by the amplitude and the 
phase variations of the relevant physical variables. This information has been 
available only piecemeal in the literature. Physical interpretations of the 
interrelationships among the detailed wave structures are given, which reveal 
important aspects of the mechanisms responsible for the energy transfer from 
the basic state to the growing waves. 


This report complements Eady’s original work, adding detailed informa- 
tion concerning the quasi-geostrophic baroclinic wave structures. 


FORMULATION AND EADY'S SOLUTION 


The derivation of the quasi-geostrophic set of equations has been amply 
documented [3] and will not be repeated here. We shall recapitulate only the 
highlights of the basic assumptions embedded in the quasi-geostrophic set of 
equations. 

Let the cartesian coordinate axes x, y, z refer to eastward, northward, 
and vertically upward directions, respectively, while u, v, w are the corres- 
ponding velocity components. A parallel flow in the x-direction of a vertically 
stably stratified Boussinesq fluid is confined within two horizontal rigid 
boundaries at z = 0 and z = H . The vertical stratification of the fluid is 
described by (the horizontally averaged) Brunt -Vais ala frequency N(z) . We 
choose U to be a characteristic horizontal velocity, L a horizontal length 
scale, L/U the time scale, and f/2 the rotation rate parallel to the z-axis. 

The geostrophic balance is achieved when the Rossby number R is 
small, i. e. , 


R q = U/(fL) « 1 


( 1 ) 


The comparable importance of stratification and rotation in the dynamical 
processes is expressed by the condition that the Froude number F(z) is of 
order unity, i.e. , 


F(z) = (f 2 L 2 )/(N 2 H 2 ) ~ 0(1) 


( 2 ) 


The hydrostatic approximation is also necessary to the derivation of the 
quasi-geostrophic equations. Note that by invoking the hydrostatic approxima- 
tion we do not allow the vertical acceleration term dw/dt in the vertical 
momentum equation; however, it is to be realized that the vertical velocity w 
itself is carried on in the formulation. Typically, the vertical velocity w is 
quite small, of the order of 6 R^U , where 6 = H/L , but because of vortex -tube 


2 


stretching the vertical velocity plays a principal role in the quasi-geostrophic 
formulations. 

Systematic scaling of the equations of motion indicates that the pressure 
deviation from the horizontally averaged hydrostatic pressure, nondimensional- 
ized by p Q UfL where is an average density for the Boussinesq fluid, is a 

nondimensional streamfunction ip for the horizontal velocities. Then the non- 
dimensional velocity components (u, v, w) are 


u = - * 


V = ip 


( 3 ) 


„ / d ,8 . 0 \ . 

w - - F -' f y K + l W K 


in which (u, v) are scaled by U, w by 5R U, (x, y) by L , z by H , and 
t by L/U. ° 


Under the above assumptions, we can write a single equation which 
expresses the conservation of potential vorticity 


( ~ 4> TT“ + $ TT - ^ h + ip + (Tip ) 1 = 0 . (4) 

y8t Y y 8x 8y ) |_/xx yy ' r z ; zj v ' 

Now, let ip= ip (y t z) be a zonally uniform, time-independent solution 
of Eq. (4) , associated with a zonal current u= -ip Q • This solution is per- 
turbed by ip* , and normal -mode solutions are assumed for the disturbances, 


and 


I * * I « I I 

Consequently, the eigenvalue problem posed for the quasi-geostrophic 
normal-mode perturbation pressure can be written as 

(u-C) [$ yy + (F? z ) z -k 2 ?] +[/3- Uyy -(Fu z ) z ]$ = 0 (6) 


with two rigid-lid boundary conditions 

(u-C)$-u$=0 on z = 0, 1 . (7) 

z z 


In Eq. (6) the dimensionless zonal wavenumber k is considered real, and the 

eigenvalue C[=C + i C ] and the eigenfunction are generally complex, 

R I 

while £ = d(R o _1 )/dy = const, represents the south-north variation of f. 

The Eady model specializes to the case in which u = z, /3 = 0, F(z) = 
const. ; Thus Eqs. (6) and (7) reduce to 

(z - C) [$ + F0 - k 2 0] = 0 . (8) 

' ' yy zz r ' ' 


and 


( z -c) $ z = 0 


( 9 ) 


/\ 1 V 

Separating out the y dependence as 0 = $ ( z) e , where I is required to be 
real, Eady obtained the closed form solutions: 
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<i>(z) = A exp(-Kz) - ( ^ ^ K C~ ) exp ( Kz ) 


( 10 ) 


* 


where 


K 



(k 2 + i 2 ) 7/3 


( 11 ) 


and the eigenvalue C is given as 


C 


-± — [V — - 

2 2K L\ 2 


COth f)(l" tanh 


ff 


( 12 ) 


where A is an arbitrary constant. 

Equation ( 12) shows that there is a critical wavelength below which all 
waves are neutral (C = 0) . The critical neutral wavenumber K is computed 
to be 1 C 


_l/ 1/ 

K c = F /2 (k 2 +i 2 ) /2 = 2.3994 . (13) 


The most unstable wave, which is associated with the maximum growth 
rate k C , occurs when 

_l/ 1/ 

F /2 (k 2 +i 2 ) /2 = 1.6062 . (14) 

The growth rate k is maximized when i = 0 , i. e. , when the disturbances 
are of infinite extent in the south-north direction. 

Note that in Eady’s model, the zonal wavenumber k and the latitudinal 
wavenumber i always appear simultaneously in the combined form (k 2 + i 2 ) . 
As a consequence, apart from the fact that i = 0 corresponds to the most 
rapidly growing mode, setting i = 0 from the outset does not cause any loss of 
generality in the dynamics in Eady’s model. However, it is worth pointing out 
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here that in attempting any improvement of Eady’s solution by incorporating 
terms of higher order in R q , neglecting the y-variation of solutions would give 

rise to inconsistency and erroneous interpretations [4], 

Employing the following shorthand notations 

= i-kc r , 

= KCj , 

= 1 + kc r * 

G^z) = (Qg 2 + Qg 2 ) exp(-Kz) -(Q^g-Qg 2 ) exp(Kz) , 

G 2 (z) = 2 Q 2 exp(Kz) , (15) 


Qi 


Qr 


Q, 


4>(z) may be rewritten as 


$(z) 




(Gj 2 + G 2 2 ) 1/z | exp [i tan” 1 (Gg/G^J 


(16) 


Once the analytical expression for the pressure perturbation 4>(z) is 
found, other dynamical variables may be derived from the quasi -geo strophic set. 
In view of the relation linking the latitudinal velocity component. 


v = R (V(.z) exp[ik(x - Ct) ]} 

© 

to the pressure perturbation, we can write 


V(z) = ik$(z) 
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The vertical component of the relative vorticity, 
& = R e (Z(z) exp[ik(x - Ct) ]} , 

may be derived from £, = V 2 ip ; thus 


Z(z) 



(Gj 


,+ a l 


'4 


•({* 


k 2 exp|i| tan“ l (G 2 /G 1 ) + 


h ! 


(18) 


Based on the thermodynamic relation in the quasi-geostrophic set, the 
dimensionless temperature perturbation, 


9 = Rg(®(z) exp[i k(x - Ct)] } , 

is given as dip'/dz. Hence, © is expressed as 


© (z) = 


( Q 2*03 2 ) 


12^/2 


(G’l + G’ 2 2 ) 


exp [ i tan -1 ( /G^ ) ] , (19) 


where 


dG dG 

G i~ nr - g 2 - nr 


Note that under the Boussinesq fluid assumption the nondimensional density per- 
turbation p is simply given as p= -0, 



We now proceed to find an expression for the vertical velocity. 


w = R e (W(z) exp[i k(x - Ct) ]} 

From the quasi-geostrophic relations, we can show that 
F" 1 W = -ik(z -C) 0 + V . 

By use of equations ( 17) and ( 19) , and rearranging, we have 


W(z) = 


qJTqJ 


) 


<■>1 + J 2> 


x k F exp[i tan" 1 (- J /J ) 

1 ^ 


where 


J 1 ' ( Z - C R )G i +C I G 2- G l • 

J 2 ’ < Z - C r) G 2- C I G 1- G 2 * 


Finally, we consider the horizontal divergence, 

div V = R (D(z) exp[i k(x - Ct) ]} 

H 6 

After some algebraic manipulation, we can show that 


D(z) = 


( Q 2 ^ Q 3 ! ) 


(M^ + M 2 2 ) 72 


xk 3 exp[i tan” 1 (- M /M ) 
J. ^ 


where 


M 


1 " < Z - V G 1 + S G 2 • 


M 2 ' < Z - V °2 - °I G 1 ' 


( 20 ) 


1 » (21) 


] • ( 22 ) 
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In summary, based on the pressure perturbation <i> in equation ( 10) , 
we have found expressions for the relevant dynamical variables in the Eady 
model. The variables are the latitudinal velocity component v, the vertical 
component of the relative vorticity t , , the temperature perturbation 0 , the 
vertical velocity w , and the horizontal divergence divV inequations (15) 
through (22). 

RESULTS AND DISCUSSION 


Firstly, the eigenvalue C(= C_ + i C ) of the Eady problem is computed. 

ft I 

Figure 1 shows the propagation speed C as a function of the scaled horizontal 

R 

wavenumber K. As was stated earlier, i = 0 , which is the most rapidly growing 

1 / 

mode, leads to the ready identification of the zonal wavenumber k= F /2 K of 

Figure 1. When the wavenumber is smaller than the cutoff wavenumber K = 

c 

2.399, C is 0.5 , implying that the unstable waves (C =£ 0, K < K ) in 
R I 

Eady’s model propagate with the same speed (which is the basic current speed 
at the mid-level z- 0.5), regardless of the wavelength. In view of the dis- 
persive nature of the real atmospheric baroclinic unstable waves, it appears 
that some of the hitherto ignored mechanisms have to be included in the Eady 
model to bring it closer to realistic atmospheric situations. 



K 

Figure 1. Propagation speed C! versus scaled wavenumber K. 

R 
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If K > K c , equation (12) shows that Cj= 0 and there are two branches 

of C , pointing to the existence of two neutral waves propagating with different 
R 

speeds, one greater than and the other smaller than the wind speed at the mid- 
level. This point will be discussed further later. 

The growth rate (K Cj) is plotted in Figure 2. The unstable region 
(K<K =2.40) and the cutoff wavenumber K beyond which neutral waves 

c o 

exist are clearly shown. The most unstable wave, which is associated with the 
maximum growth rate, is seen to have wavenumber K = 1.6062 . It is also 
apparent that in the unstable region (C J ^ 0 , K < K c ) , for every growing mode 

solution corresponding to Cj> 0 there exists a damped wave with Cj < 0 . 



Figure 2. Growth Rate KCj versus scaled wavenumber K. 

We now present detailed information about the wave structure for the 
disturbances of varying zonal wavelengths. Since an eigenfunction is generally 
complex, the amplitude (or the modulus) and the phase will be displayed 
separately. For example, in the case of the perturbation pressure $ (z) , 
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1 / 

the amplitude, [A/(Q 2 + Q 2 ) ] [(G 2 + G 2 ) z ] , and the negative value of the 

phase, -tan _1 (G /G ), in accord with equations (15) and (16) are shown in 

Figures 3(a) and 3(b), respectively. Note that the eigenfunction solution of 
equation (8) contains an arbitrary constant A , which is set to be unity; there- 
fore, the amplitude plots display only the relative magnitude of variation with z . 
Also, for the sake of better graphical presentation, different scale factors (for 
the plot of amplitude variations) and phase shift factors (for the plot of phase 
variations) are employed for each variable in the figures. 



(a) Amplitude variations. Scale factors are 0.5 for | <f> | , 0.5 for 
| ©| , 1.0 for |W| , and 0.3 for |D| . 



(b) Negative relative phases of perturbation pressure <i> and 
temperatime ©. Phase shift factor is 11.35°. 



Z=1 


z=o 


(c) Negative relative phases of vertical velocity W and 
divergence D. Phase shift factor is 11. 35° . 


Figure 3. Amplitude and phase variations of the most unstable mode 

(K = 1.6061, C = 0.5, C = 0.1929). 

£\ I 
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The wave structures of the most unstable mode are described in Figure 
3. These plots may also be construed as being qualitatively representative of 
the wave structures of unstable waves in general. Because of the simple linear 
wind profile and the rigid boundary conditions, the amplitudes of the eigen- 
functions of the Eady unstable waves are symmetrical with respect to the mid- 
level. The temperature and pressure perturbations are largest at the boundaries. 
The vertical velocity | W(z) | is zero at both boundaries and reaches the maxi- 
mum at the mid-level. On the other hand, the horizontal divergence | D( z) | is 
a maximum at both boundaries and decreases to its minimum at the mid-level. 
Figures 3(b) and 3(c) show the vertical structure of the negative phase plot of 
the dynamical variables of the most unstable mode. The negative phase dia- 
grams show the actual side view of the wave structures. Combining the ampli- 
tude and the negative phase information, the x-z cross section over one 
wavelength is constructed in Figure 4 to illustrate the interrelations among 
dynamical variables. In Figure 4, the size of the letterings indicates the rela- 
tive magnitude of the variables denoted. Both the pressure trough (high 
pressure, denoted by 4> + ) and the pressure ridge (low pressure, denoted by ) 

are tilted upward and backward. The warm tongue (the perturbation temperature 
maximum, denoted by © ) lies ahead of the pressure trough but to the rear of 

the pressure ridge. Again, the cold tongue (the perturbation temperature mini- 
mum, denoted by © ) is located at half-wavelength apart from the warm tongue. 

In a similar fashion, the vertical velocity distribution W, denoted by vertical 
arrows, and the horizontal divergence D, denoted by horizontal arrows, are marked 
in Figure 4. At the same time, the south-north (y) component of velocity V 


-© 




Q + * + 


O -z»i 


0 + ® 


©I e_ $>- ®fe+ <i*. 

©_J© 4 >_ 

e_ I© $>_ e+ 1® 

— 9 _ O =^ 4 >_ 0 + © 


$ + 




4> + -z-o 


Figure 4. The x-z plot showing the interrelationships among the dynamical 
variables. The size of the letterings indicates the relative 
magnitude of the variables. 
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expressed in equation ( 17) , which is shown by (x) ,0 in Figure 4, has the same 
amplitude variation as the perturbation pressure <3? , but the phase of V is 
7r/2 ahead of the phase . The x-z cross section of Figure 4 reveals the well- 
known potential energy release mechanism associated with the baroclinically 
unstable waves [5]. If the kinetic energy of the wave field is to increase with 
time, the vertical velocity has to be positively correlated with the perturbation 
temperature so that on the average warm fluid rises while cold fluid falls. The 
correlation between w and 6 in Figure 4 for unstable waves verifies the 
preceding condition. Another significant aspect of the baroclinically unstable 
wave energetics is that the correlation between the south-north velocity com- 
ponent v and the perturbation temperature 6 should be positive, suggesting 
that the unstable waves, on the average, transport heat down the basic tempera- 
ture gradient in the south -north direction. This is also shown in Figure 4; 
warm (cold) fluid is correlated with the velocity component into (out of) the 
plane of paper. In summary, the detailed wave structures of Figure 3, together , 
with the x-z plane plot, show how baroclinically unstable waves withdraw energy 
from the store of available potential energy of the basic current which is sus- 
tained by the latitudinally varying temperature field. 

Now, it will be of interest to examine in detail the peculiar features of 

the baroclinic wave structures as the wavenumber K increases. Figure 5 

illustrates the structure of an unstable mode at wavenumber K = 2.2 , which has 

a reduced growth rate (C = 0.5, C = 0.0961) . Although the overall 

R I 

patterns of the wave structure at this wavenumber are quite similar to those of 
the most unstable mode, some differences are easily discernible with respect 
to the negative phase diagrams of Figures 5(b) and 5(c). Whereas the phase 
line of © is further stretched toward the horizontal, the phase line of <£ is more 
erect. Perhaps, the behavior of the vertical velocity W is the most interesting 
to analyze because other variables may be readily expressed in terms of W in 
the quasi-geostrophic formulation. The constant phase line of W of Figure 
5(c) still slopes upward and backward, which is characteristic of unstable 
baroclinic waves; however, the slops is steeper (toward the vertical) than the 
case of the most unstable mode of Figure 3(c) . This steepening of the phase 
line of W will become more apparent as the wavenumber approaches the 
cutoff wavenumber K . 

Figure 6 shows the structure of an unstable wave with the wavenumber 
closer to the cutoff wavenumber. At this wavenumber (K= 2.35) , the phase 
speed is still 0. 5, but the growth rate is quite small (Cj = 0.0476) . The phase 

lines of © and D are further stretched toward the horizontal, and a trend may 
be recognized in that the phase lines of © and D tend to a step-like shape. 
However, the phase lines of <3? and W are further steepened toward the vertical. 
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(c) Negative relative phases of vertical velocity W and 
divergence D. Phase shift factor is 5. 13°. 

Figure 5. Amplitude and phase variations of an unstable wave (K= 2.2, 

C = 0. 5, C = 0.0961). The wavenumber K(= 2.2) is greater than 
R I 

the wavenumber of the most unstable mode K= 1.6061. 
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1 0 1 

(a) Amplitude variations. Scale factors are 0. 5 for | $ | , 0.3 for 
| ©|, 0.3 for | W|, and 0.1 for |D|. 



Z=1 


0 


z=o 

It 


(b) Negative relative phases of perturbation pressure $ and 
temperature ©. Phase shift factor is 2.49°. 



(c) Negative relative phases of vertical velocity W and 
divergence D. Phase shift factor is 2.49°. 


Figure 6. Amplitude and phase variations of an unstable wave (K= 2.35, 
C^= 0.5, Cj= 0.0476). The wavenumber K(= 2.35) is close 
to the critical wavenumber K_. 
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When the wavenumber equals the critical wavenumber, we have a 

marginally stable mode. The structure of this marginally stable mode (K = 

2.3993, C = 0.5, C = 0) , as shown in Figure 7, exhibits features which are 
R I 

the limiting behavior of the unstable waves. For this marginally stable mode, 
the amplitude functions of Figure 7(a) are still symmetric with respect to the 
mid-level. Both the temperature perturbation © and the horizontal divergence 
D vanish at the mid-level. The phase lines of $ and W are now vertical, and 
the phase lines of © and D are step functions with steps at the mid-level. These 
vertical phase line patterns are typical of the neutral waves (C^= 0) in general, 

where the mechanism for converting energy from the basic current to perturba- 
tions is absent. Tokioka [6], by computing the energy transfer rate in terms 
of the vertical velocity, demonstrated that the phase line of W is vertical [i. e. , 
da /dz = 0, if W is written as W = |W| exp(ia ) ] when there is no energy release 
from the basic state to the disturbances. 

In the stable wave region where the wavenumber is larger than the cutoff 
wavenumber, as is shown in the eigenvalue plots of Figures 1 and 2, there 
appear two neutral modes (C^. = 0) , each of which has a different phase speed 

C , one larger than and the other smaller than the basic current speed at the 
R 

mid-level. 

The wave structure of a fast-propagating (C_. > 0.5) neutral mode (K = 

R 

3.5, C = 0.7131, C = 0) is shown in Figure 8. For this neutral mode, th,e 

amplitude functions of Figure 8(a) are no longer symmetric with respect to the 
mid-level. It is also evident that the amplitude distribution functions are 
enhanced in the region near the upper boundary. The perturbation temperature 
© is zero at a level in the lower half of the fluid depth, and the phase line of © 
has a step at this level. The horizontal divergence D becomes zero at a level 
between the mid-level and the upper boundary, and the phase line of D has a 
step at this level. As was pointed out earlier, the phase lines of $ and W are 
vertical, allowing no energy conversion mechanism from the basic state to the 
disturbance. 

The wave structure of a slow-propagating (C^ < 0.5) neutral mode (K = 

R 

3. 5, C = 0. 2869, C = 0) is shown in Figure 9. The bulk of the disturbance is 
R I 

confined to the bottom half -depth of the fluid. This slow-propagating neutral 
mode exhibits features opposite to the fast -propagating neutral mode described 
in Figure 8. 
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(a) Amplitude variations. Scale factors are 0. 5 for | $ | , 0.2 for 
| ©|, 0.25 for | W|, and 0.1 for |D|. 



( b) Negative relative phases of vertical velocity W and 
perturbation pressure . Phase shift factor is 0. 0. 



(c) Negative relative phase of perturbation temperature ©. 
Phase shift factor is 0.0. 



(d) Negative relative phase of divergence D. 
Phase shift factor is 0.0. 


Figure 7. Amplitude and phase variations of a marginally stable mode 
(K = 2.3996, C R = 0.5, = 0.0). 
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(a) Amplitude variations. Seale factors are 0.5 for | <i> | , 0.1 for 
| ©| , 0.2 for | W | , and 0.05 for |D| . 



(b) Negative relative phases of perturbation temperature © and 
divergence D. Phase shift factor is 0.0. 



(c) Negative relative phases of vertical velocity W and perturbation 
pressure $ . Phase shift factor is 0. 0. 

Figure 8. Amplitude and phase variations of a fast-propagating stable 
mode (K = 3.5, C R = 0.7131, C = 0.0). 
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(a) Amplitude variations. Scale factors are 0.5 for | $ | , 0.1 for 
I© |, 0.2 for |W|, and 0.05 for |D|. 



(b) Negative relative phases of perturbation temperature © and 
divergence D. Phase shift factor is 0. 0. 



(c) Negative relative phases of vertical velocity W and 
perturbation pressure <£. Phase shift factor is 0.0. 


Figure 9. Amplitude and phase variations of a slow-propagating stable 
mode (K = 3.5, C = 0.2869, C = 0.0). 

JA 1 
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A physical interpretation is helpful in explaining the transition from the 

neutral modes to the unstable modes as the wavelength increases. In the stable 

region (K > K ) the horizontal length scale (i.e. , wavelength) L is small so 
C s 

that the corresponding penetration scale height [7,8], of the order of fL /N, 

s 

is smaller than the half -depth of the fluid. As a consequence, the motion forced 
by the external conditions at the rigid boundaries cannot penetrate fully into the 
interior of the fluid. The resulting motion is thus concentrated in regions 
adjacent to the boundaries. Under these circumstances, the two neutral modes, 
energetically incapable of converting potential energy from the basic current to 
the disturbances, are simply advected by the basic current. Since the basic 
current is linearly increasing with z, the mode whose motion is concentrated 
near the upper boundary is advected faster than the mode whose motion is 
concentrated near the bottom boundary. 

As the wavelength increases to the point where the penetration scale 
height becomes comparable to the half-depth of the fluid, the two neutral modes 
coalesce to form a single mode, and instability begins to set in. 
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